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Abstract. We discuss how orbital degeneracy, which is usually removed by a 
cooperative Jahn- Teller distortion, could under appropriate circumstances lead rather 
to a Resonating Valence Bond spin-orbital liquid. The key points are: i) The tendency 
to form spin-orbital dimers, a tendency already identified in several cases; ii) The 
mapping onto Quantum Dimer Models, which have been shown to possess Resonating 
Valence Bond phases on the triangular lattice. How this program can be implemented 
is explained in some details starting from a microscopic model of LLN1O2. 



PACS numbers: 75.10.Jm, 05.50.+q, 05.30.-d 
1. Introduction 

It is very useful to classify quantum magnets according to the symmetry (if any) that 
is broken in the ground state. When the SU(2) symmetry is broken, the system 
usually sustains some kind of long-range magnetic order, although some more exotic 
examples involving quadrupolar order have been recently discussed P] E]. When the 
SU(2) symmetry is not broken, a translation symmetry may or may not be broken 
depending on the lattice topology. The standard case in which no lattice symmetry is 
broken is that of systems in which it is possible to build a singlet inside the unit cell, 
as for instance in a ladder 0. If however this is not possible, as in all systems with 
half-integer spins and an odd number of sites per unit cell, the simplest way to keep 
SU(2) symmetry is to break the translation symmetry so that the new unit cell contains 
an even number of sites per unit cell. The typical example is the S=l/2 J\ — J2 chain, 
which has been explicitly shown quite some time ago by Majumdar and Ghosh [4] to 
have two-fold degenerate dimerized ground state when J 2 = Ji/2. 

Another possibility to keep SU(2) symmetry without breaking any lattice symmetry 
has been put forward by Anderson[5j in 1973. Concentrating on S=l/2 magnets, he 
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suggested that under appropriate circumstances the ground state might be a linear 
combination of valence bond (VB) states, i.e. states in which sites are paired to 
form singlets. Clearly each individual state breaks at least some of the translational 
symmetries of the underlying lattice, but the translational symmetry is restored by 
the superposition of such valence bond states. Such a wave function is known as a 
Resonating Valence Bond (RVB) state. 

To identify such a ground state in realistic models of quantum antiferromagnets 
turned out to be much more difficult than anticipated. The prediction of Fazekas 
and Anderson]?)] that this might be the case for the S=l/2 Heisenberg model on the 
triangular lattice, based on estimates of the energy of ordered and valence bond states 
including perturbation corrections, is not supported by recent numerical investigations, 
which all point to 3-sublattice antiferromagnetic order [7 . 

The most serious candidate still around is the S=l/2 Heisenberg antiferromagnet 
on the kagome lattice. No evidence of magnetic long-range order could be found so 
far, and the proliferation of low-lying singlets observed in exact diagonalizations of 
finite clusters jH] can be fairly well described in the short-range RVB subspace of valence 
bond states involving only nearest neighbours 0. However, several treatments based on 
some effective Hamiltonian have reached the conclusion that the ground state support 
some kind of valence bond order ^UJ El E] and hence breaks translational symmetry. 
Unfortunately, the resulting unit cell is so large, and the singlet-singlet gap accordingly 
so small, that this prediction cannot be cross-checked by the only unbiased numerical 
approach available so far, namely exact diagonalization, and the issue is likely to remain 
open for quite some time. 

In fact, it has only been possible so far to unambiguously identify an RVB ground 
state in a very minimal description of fluctuations in the RVB subspace that goes 
under the name of Quantum Dimer Model (QDM)^3|. In this approach, valence bond 
configurations are assumed to build an orthogonal basis of the Hilbert space, and the 
effective Hamilonian contains kinetic terms that shift dimers around loops and potential 
terms that favour or penalize specific local configurations of dimers. For the triangular 
lattice, the simplest model is defined by the Hamiltonian: 

H = vJ2(\-)iS\ + \//)(//\) -*£ ( I J7 > < Z;/ I + I /:/ > < _T I ) 

where the sum runs over all plaquettes including the three possible orientations. The 
kinetic term controlled by the amplitude t changes the dimer covering of every flippable 
plaquette, i.e., of every plaquette containing two dimers facing each other, while the 
potential term controlled by the interaction v describes a repulsion (y > 0) or an 
attraction (v < 0) between dimers facing each other. Since a positive v favors 
configurations without flippable plaquettes while a negative v favors configurations with 
the largest possible number of flippable plaquettes, one might expect a phase transition 
between two phases as a function of v/t. The actual situation is far richer though. 
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As shown by Moessner and Sondhi[T3], who calculated the temperature dependence of 
the structure factor, there are four different phases: i) A staggered phase for v/t > 1, 
in which the ground-state manifold consists of all non-flippable configurations; ii) A 
columnar ordered phase for v/t sufficiently negative; iii) An ordered phase adjacent 
to it which probably consists of resonating plaquettes which make a 12-site unit-cell 
pattern [T5]: iv) A liquid phase with a featureless and temperature independent structure 
factor. This last phase has been interpreted as a short-range (RVB) phase in which all 
correlations decay exponentially at zero temperature, an interpretation confirmed by 
recent Green's function Quantum Monte Carlo simulations|16j. which have established 
the presence of topological degeneracy in this parameter range, a clear characteristic of 
the RVB phase. 

This result defines a new line of research in the field: Indeed, rather than 
investigating directly the properties of the Heisenberg model on a given lattice, one 
can try to identify models for which a VB subspace is a reasonable variational subspace, 
derive an effective QDM, and study it along the same lines as the minimal model on 
the triangular lattice. In that respect, a natural candidate is the trimerized spin- 1/2 
Heisenberg model on the kagome lattice. An effective model in terms of the total spin a 
and a chirality pseudo-spin r per strong triangle has been derived fT7\ ITS]. It is defined 
on the triangular lattice built by strong triangles and can be written [T9]: 

K = t £ ** • - ^ • W ~ 4e ^ • f i)> W 

where J' is the weak coupling between the strong triangles of the trimerized lattice, 
and where the vectors eV,- have to be chosen among e\ = (1,0), e 2 = (— | , — 2 )' 
e*3 = (— |, ^) according to the pattern of Fig. [T] A mean-field decoupling of spin and 
chirality has identified nearest-neighbour valence bond states as the lowest solutions j!8j. 
and, following Rokhsar and Kivelson^SI, a QDM has been derivedj20]- Unfortunately, 
the properties of this model could not be studied so far. First of all, they involve kinetic 
terms on longer loops than the above-mentioned minimal model, but more importantly, 
it is impossible to formulate it in such a way that all off-diagonal matrix elements 
are negative, so that the Quantum Monte Carlo methods that were successful for the 
minimal model cannot be used. 

The effective spin-chirality model has another remarkable feature though: It is 
formally very similar to the spin-orbital models that are used to describe Mott insulators 
with orbital degeneracy. Indeed, as discussed in great details by Kugel and Khomskii[2*T|. 
when the local symmetry is such that different orbitals can be occupied in the open 
shell of the magnetic ions of a Mott insulator, this extra degree of freedom can be 
described by a pseudo spin, and the resulting model is roughly speaking of the same 
form. Given the very different precise forms this model can take for specific systems, a 
general discussion cannot be attempted here. Rather, we concentrate in the next section 
on a specific example of Mott insulator with orbital degeneracy for which we believe 
that RVB physics might be realized. More general comments will be given in the last 
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Figure 1. Triangular lattice on which the spin-chirality Hamiltonian is defined. The 
unitary vector for the bond is indicated by solid lines (e) t = ei), dashed lines [e^ = e^), 
and dotted lines (e M = e^). 

two sections of the manuscript. 

2. The spin-orbital model of LiNi02 

The Mott insulator LiNi02 and its cousin NaNiC>2 are isostructural and isoelectronic. 
The crystal structure can be envisaged as a sequence of slabs of edge sharing octahedra 
of oxygen O 2- ions. Metal ions sit at the centers of octahedra. There are two kinds of 
slabs: in A slabs, at every center of octahedra there is a Ni 3+ , whereas in the B slabs, 
one finds either Li + or Na + ions. A and B slabs alternate (see Fig. |2J). The Ni ions 
form well-separated triangular planes. The Ni 3+ ions are in the S = 1/2 low-spin state, 
which allows for twofold orbital degeneracy between the d 3z 2_ r 2 and d x 2_ y 2 orbitals (see 
Fig. ED. 




Figure 2. Left: AMO2 structure. Ni ions are located in the middle of the O 
octahedra. Right: Local structure and degenerate orbitals: d 3z 2_ r 2 (left octahedron) 
and d x 2_ y 2 (right octahedron) 

Surprisingly enough, the two systems have very different properties: NaNiC>2 
undergoes a high temperature Jahn- Teller distortion, followed at low temperature by 
the antiferromagnetic ordering of ferromagnetic triangular planes, which makes it a 
standard example of orbital ordering followed by magnetic ordering|22|- By contrast, 
no ordering could be detected in LiNi02, and some kind of freezing seems to take place 
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below 8 Kj2H]- The best samples of LiNi02 are non stoichiometric however, some Li sites 
being occupied by Ni atoms, and this has been invoked to explain the difference. The 
trend upon approaching stoichiometry is not clear though, and this striking difference 
between the two systems, in particular the lack of any sign of a cooperative Jahn- Teller 
transition in LiNi02, suggests to look for alternative explanations. In the following, 
starting from a realistic microscopic description of the system, we study the possibility 
to stabilize an RVB spin-orbital liquid in the absence of any disorder. 



2.1. Microscopic Model 



A fairly general description of this system is given by a Kugel-Khomskii Hamiltonian 
defined in terms of Wannier functions centered on the Ni sites by two hopping integrals 
t h and t' h , the on-site Coulomb repulsion U and the Hund's coupling J which, on a given 
bond, takes the form 24j 
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with the usual definitions for the projectors on the singlet and triplet states of a pair of 
spins: 
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The vectors nf- depend on the orientation of the bonds and are given by n\ 2 = (0, 0, 1), 

(— ^r,0, — |) for the 3 orientations respectively. The 
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operators Tj are pseudo-spin operators acting on the orbitals. For the local geometry 
shown in the right panel of Fig. El th and t' h correspond to the hopping between pairs 
of G? 32 2_ r 2 and d x 2_ y 2 respectively, the hopping between the d% z 2_ r 2 on one site and the 
d x 2_ y 2 on the other being zero by symmetry. Of course, all bonds are equivalent, but once 
a basis, i.e. a pair of local orbitals, has been chosen, the Hamiltonian takes a different 
form on bonds with different orientations. Note that other forms of the microscopic 
Hamiltonian including explicitly O orbitals have been used|251 123 HZj with somewhat 
different conclusions. 



2.2. Mean-field phase diagram 

Inspired by the results obtained on the trimerized kagome modelJTH], spin and orbital 
degrees of freedom can be decoupled in a mean-field way [21] ■ This leads to a phase 
diagram in which phases can be distinguished by the mean value of the orbital and/or 
spin part of the Hamiltonian on each bond. The resulting phase diagram is remarkably 
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Figure 3. Mean-field phase diagram on a 16-site cluster as a function of hopping 
integral versus Hund's coupling. The grey phase is the ferromagnetic phase, with the 
classical phase boundaries between different types of orbital ordering. 

While the planes of NaNiC>2 are known by now to be ferromagnetic, suggesting 
that NaNiC>2 is in one of the ferromagnetic phases, LiNi02 is expected to be in one of 
the antiferromagnetic phases. The orbital and spin structure of the antiferromagnetic 
phases is depicted in Fig. HI except phase A, which consists of SU(4) plaquettes[28 
and cannot be described along these lines. Since niether orbital nor magnetic ordering 




D E 

Figure 4. Spin and orbital structure in the singlet phases of the mean-field phase 
diagram. Solid line indicates AF, dashed line FM spin correlations. 

could be detected in LiNi02, and since phases B and D have a simple orbital ordering 
pattern while phase E is likely to be antiferromagnetically ordered, let us concentrate on 
phases C and C Both phases are characterized by strong dimer singlets forming different 
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regular dimer coverings of the triangular lattice. On each dimer the orbitals are parallel, 
and they correspond to d 3z 2_ r 2, d 3x 2_ r 2 or d 3y 2_ r 2 depending on the orientation of the 
bond. Note that all these orbitals are Jahn- Teller active, leading in all cases to two 
long and four short Ni-0 bonds. One might be tempted to conclude that these phases 
correspond to two types of valence bond solids with the patterns depicted in Fig. |U 
The mean-field approach has a very remarkable property however: In addition to the 
mean-field solutions with lowest energy shown in Fig. HJ there are several other mean- 
field solutions of the self-consistent equations with energies very close to the lowest 
energy corresponding to other dimer coverings of the triangular lattice In such 
circumstances, going beyond mean-field is likely to couple these solutions, and the 
relevant model would then be a QDM describing resonances between these states, a point 
of view favoured by exact diagonalizations of finite clusters. So it is more appropriate 
to think of these phases as a region of parameters where all dimer coverings are relevant 
states for low-energy physics. 

2.3. Effective Quantum Dimer Model 

Starting from all dimer configurations mentioned in the previous section, a QDM has 
been derived jHHj- It involves a competition between kinetic processes and dimer-dimer 
repulsion. A miminal version of the model is defined by: 

H= - t?,(\//){-\ + h.c))- t'E^ 

+ ^e(|//)(//I + I_-)(_ _ |), 

where the sums run over the 4-site and 6-site loops with all possible orientations. 
Although the repulsion is a higher order process, hence quite small, and although 
the ratio t'/t is in principle fixed by the perturbative expansion, these parameters are 
treated as free to make contact with the Rokhsar-Kivelson model on the triangular 
lattice. The main difference with the effective model derived for the trimerized kagome 
antiferromagnet is that the off-diagonal elements are now all negative. This is in practice 
extremely important since it allows one to use Quantum Monte Carlo simulations. 

The phase diagram of the model has been derived using exact diagonalizations of 
finite clusters and Green's function Quantum Monte Carlo [30J. As shown before, the 
most convenient way to identify an RVB phase is to look for topological degeneracy since 
it is at least partially lifted in all other phases. The resulting phase diagram is shown 
in Fig. E| It contains a large RVB liquid phase which connects the relevant parameter 
range for LiNi02 (t'/t ~ 2 and V small) to the RVB liquid phase of the minimal model 
(t' = 0). Translated into spin-orbital language, this RVB phase corresponds to a spin- 
orbital liquid wiht no symmetry breaking and no phase transition, in agreement with 
the phenomenology of LiNi02- 




+ h.c. 



(4) 
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Figure 5. (Color online) Phase diagram in the t'—V plane. A wide disordered region 
extends all the way from the standard QDM it' ft = axis) to the purely kinetic QDM 
(V/t — axis). The description of the symbols is given in the text. 



3. Discussion 

Let us now comment on how generic the mechanism proposed in the context of the 
spin-orbital model of LiNi02 might be. The main ingredients to get an RVB spin- 
orbital liquid are: 1) The spontaneous formation of dimers; 2) The degeneracy or 
quasi-degeneracy of the energies of the wave-functions constructed out of different dimer 
coverings; 3) The presence of an RVB phase in the relevant QDM. Let us comment on 
these points separately. 

The tendency of spin-orbital models to spontaneously form dimers is well 
documented. The possibility to stabilize dimerized ground states due to orbital 
degeneracy has first been put forward by Feiner et al in the context of a realistic 3D 
model|2T|. Shortly after, the presence of a dimerized ground state has been explicitly 
proven for a simple minimal ID model by Kolezhuk and Mikeska|22]; a result generalized 
shortly after by Pati et a/j33j in the context of a model defined by the Hamiltonian: 



# =£ [Ji&Jj + J 2 T i .f j + K(§ i J j )(f i .'r j j\ (5) 

(iJ) 

with K > 0. When J±/K — J2/K = 3/4, it can be rewritten as 

H = KJ2&.S 3 + 3/4) + 3/4) (6) 

(id) 

Each term is obviously positive, and since Si.Sj + 3/4 (resp. Tj.Tj + 3/4) is the projector 
on the spin (resp. orbital) triplet, the two wave-functions with alternating spin and 
orbital singlets are zero energy eigenstates, hence ground-states [32] • Pati et al have 
shown that this dimerized phase extends to a very large portion of the phase diagram 
around this point. From that point of view, the identification of dimer phases in the 
context of spin-orbital models of BaVS3[29j and of LiNiC^jll] is not unexpected, and 
the tendency to dimerize can be considered to be a generic trend of spin-orbital models. 

What seems to be more specific to these spin-orbital models of LiNi0 2 and BaVS 3 
is the quasi-degeneracy of all nearest-neighbour dimer coverings. But in fact, this can 
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be traced back to a rather generic feature of spin-orbital models, namely the remarkable 
symmetry properties of the orbital part of the Hamiltonian. As can be clearly seen 
in the spin-orbital model of LiNi02, the orbital part does not have the same form in 
the three directions of the triangular lattice, a property encoded in the n?j vectors. 
So the Hamiltonian is only invariant if one simultaneously performs the same rotation 
in real space and in pseudo-spin space. For purely orbital models, this is known to 
have remarkable consequences |34[ 1351 For spin-orbital models, this implies that 
dimers with different orientations involve different orbital wave-functions, as can be 
clearly seen in phases C and C. What controls the energy of a given dimer configuration 
is then the residual dimer-dimer interaction. It turns out that simple patterns having 
all dimers parallel to each other are not naturally favoured if the anisotropy of the 
orbital part is strong because it is impossible to gain energy in the other directions. On 
a lattice such as the triangular lattice with a large connectivity, it is then much more 
favourable to adopt configurations where dimers are not parallel to each other. The 
energy difference between such configurations however is not really significant, and it is 
better to look at such states ClS db variational basis. 

Finally, too little is known at that stage about RVB phases in QDM to draw 
general conclusions, but it seems plausible that the presence of an RVB phase between 
two valence-bond phases is the generic alternative to a first order transition. 

To summarize, the tendency toward dimerization is a rather general feature of spin- 
orbital models. When confronted to a lattice such as the triangular lattice, for which 
QDM's are known by now to possess RVB phases, there is a real chance for quantum 
fluctuations to stabilize an RVB ground state. 

4. Conclusion 

We have shown that orbital degeneracy can, under special but neither unrealistic nor 
fine-tuned conditions, lead to a spin-orbital RVB ground state. Clearly, this is not the 
most common situation. Indeed, orbital degeneracy usually leads to a cooperative Jahn- 
Teller transition, resulting into an effective spin Hamiltonian with a symmetry different 
from that of the original lattice[2I]. This is not either the only route to spin-orbital 
liquid behaviour [3H]- However, the tendency of spin-orbital models to spontaneously 
dimerize is strong enough to make this a promising route towards RVB physics. Whether 
LiNi02 is the first example remains to be seen. To make progress on this issue will 
require not only further theoretical work to better understand the relevant microscopic 
model and its possible connection to a QDM, but also and maybe more importantly 
further experimental investigations to unambiguously identify the orbital and magnetic 
properties of stoichiometric samples. 

We acknowledge useful discussions with M. Ferrero and D. Ivanov. This work was 
supported by the Swiss National Fund and by MaNEP. 
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